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INTRODUCTION 


Experimental  determinations  of  the  temperature  T  dependence  of  the 
electrical  resistivity  p  in  amorphous  metals  have  generally  been  made  at 
constant  pressure  P.  Nevertheless,  most  theoretical  studies  of  electrical 
transport  in  amorphous  metals  have  treated  the  variation  of  p  with  T  at 
constant  volume  V.  This  has  been  the  case  because:  (1)  there  were  few 
published  measurements  of  p(T)  at  pressure;  (2)  the  available  data  did  not 
seem  to  show  any  regularity;  and  (3)  the  early  experimental  studies  (refs  1,2) 
which  were  performed  in  high  p,  TM-based,  amorphous  alloys,  exhibited  pressure 
coefficients  of  resistivity  that  were  considerably  smaller  in  magnitude  than 
those  in  lower  p,  crystalline  alloys  (ref  3). 

However,  recent  experimental  studies  of  the  effects  of  pressure  on 
electrical  transport  in  amorphous  metals  by  Fritsch,  Wilier,  Wildermuth,  and 
Luscher  (ref  4)  Indicate  that  thermal  expansion  effects  cannot  generally  be 
ignored.  Recently,  Hafner  (ref  5)  found  that  thermal  expansion  effects  are 
important  in  determining  the  T-dependence  of  p  at  constant  pressure  in 
a-Mg7Zn3  in  a  pseudopotential  based,  diffraction  model  calculation.  Further¬ 
more,  reexamination  of  the  early  data,  assuming  reasonable  values  for  the  bulk 
moduli  B-j,  suggests  that  thermal  expansion  effects  can  also  be  important  in 
high  resistivity  alloys. 

Thus,  it  is  our  objective  to  give  a  simple  description  of  the  way  that 
thermal  expansion  effects  influence  the  p(T)  in  amorphous  metals.  Although  we 
present  examples  in  the  context  of  the  diffraction  model,  our  principal 
results  are  not  transport  model  specific.  Our  approach  is  to  assume  that 

l 

References  are  listed  at  the  end  of  this  report. 
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macroscopic  properties  such  as  By(T),  the  Gruneisen  constant  y(T) , 

3inp/3£nV | ^ ,  etc-,  are  given  and  that  the  thermal  expansion  (i-e.,  V(T))  is 
described  by  Gruneisen' s  law  (ref  6).  This  approach  is  in  contrast  to  the 
tour  de  force  of  Hafner  (ref  5),  in  which  all  relevant  quantities  are  computed 
from  first  principles  using  pseudopotential  methods  specifically  for  a-KgjZng 
and  a-CayMg3- 

THEORY 

The  variations  of  P  with  volume  variations  produced  by  thermal  expansion 
are  small  compared  to  p  for  0  <  T  <  23  in  amorphous  metals.  Thus,  one  may 
expand  p  in  a  Taylor  series  in  in(V)  and  neglect  terms  higher  than  first 
order: 

p(T,V)  =  p(T,V0)[l  +  Q(T,V)]  (la) 

where  the  thermal  expansion  part, 

3  in  P I 

Q(T,V)  =  *n(V/V0)  - - 1  *  (lb) 

3  inV  |  T ,  VQ 

and  V0  is  the  mean  atomic  volume  at  0°K  and  one  atmosphere.  Similarly,  the 
isochoric  variations  in  p  are  small  compared  to  p  for  0  <  T  <  20  so  that  we 
may  rewrite  Eq.  (1)  as 

P(T, V)  =  Po[l  +  R(T,V0)  +  Q(T, V) ]  (2) 

where  pQ  =  p(0,VQ)  and  R(T,V0)  is  the  isochoric  relative  change  in 
resistivity,  i.e.  , 

R(T,V0)  =  p(T,V0)/p0  -  1  (2a) 

Note  that  we  drop  the  higher  order  R(T,VQ)  Q(T,V)  term  to  obtain  Eq .  (2). 

Now,  4n(V(T)/Vo)  «  V(T)/V0  -  1  since  V(T)/V0  -  1  «  1  for  0  <  T  <  29. 
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Thus,  Gruneisen's  law  (ref  6)  of  thermal  expansion  yields  at  constant 


pressure : 

£n( V/VQ)  =  y!(T)7btV0  (3) 

where  e(T)  is  the  mean  thermal  energy,  Y  is  the  Gruneisen  constant,  and  BT  is 
the  bulk  modulus.  The  mean  thermal  energy  can  then  be  evaluated  for  a  Debye 
model,  yielding 

THT  =  3NkB8F(T/9)  (4a) 

where  N  is  the  avagodro's  number,  kg  is  the  Boltzmann  constant,  and 


F(x)  =  x4/  3z3dz/(ez-l) 
0 


The  limiting  forms  of  F(T/0)  are 

F(T/0)  =  (  tt4/5)  (T/  9) 4  for  T/ 0  «  1  (4c) 

and 

F(T/0)  =  T/0  for  T/0  »  1  (4d) 

Figure  1  exhibits  graphs  of  F(T/0)  and  F’(T/0)  versus  T/0.  F(T/8)  and  F’(T/0) 

% 

describe  the  T-dependence  of  e(T)  and  the  lattice  specific  heat  CV(T), 
respectively,  in  the  Debye  model.  We  shall  see  below  that  the  thermal 
expansion  part  of  the  isobaric  p(T)  is  proportional  to  F(T/ 0)  and  that  the 
thermal  expansion  coefficient  is  proportional  to  F'(T/0). 

Thus,  our  principal  result  at  constant  pressure  is 

p(T)  =  p(T,V(T))|p  =  pQ(l  +  R(T,V0)  +  Q(T) )  (5) 

where 

Q(T)  =  Q(T,V(T) ) |p  =  K(T)F(T/0)  (6a) 

with  K(T)  =  (3NkB0(T)Y(T)/Bx(T)Vo)( 34np/3inV)  (6b) 

T,Vq 

and  R(T,VQ)  is  the  isochoric  relative  change  in  P  defined  in  Eq.  (2a).  The  T- 
dependence  of  K(T)  is  much  weaker  than  that  of  F(T/ 0) .  Moreover,  at  low  T, 
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where  the  temperature  variation  of  K(T)  is  largest,  F(T/9)  varies  very  rapidly 
with  T  (i.e.,  like  T4)  and  is  small.  Thus,  it  is  a  good  approximation  to 
neglect  the  T  variation  of  K(T)  and  to  use  K(T)  K  K( 9) . 

An  alternate  expression  for  K( 9)  can  be  derived  in  terms  of  the  thermal 
expansion  coefficient  0V(T), 

BV(T)  5  ( 3 &nV / 3T ) p  -  3NkBy(T)F' (T/ 9)/BT(T)VQ  (7) 

where  we  have  neglected  the  e(T)  3(  y/BxV0) /  ST  term  and  a  term  arising  from 

30/3T.  (Equation  (7)  is  equivalent  to  Eq .  (9)  of  Reference  5.)  Then 

substituting  Eq .  (7)  in  Eq.  (6b)  at  T  =  9,  one  finds 

K(8)  =  0  ( 9) 9( 9) ( 3£np/ 3&nV)  /F’(l) 

e>Vo 

-  1.05  90v( 0)( 3Znp)/ 3£nV)  (8) 

Mo 

(Equivalent  expressions  can  be  given  if  0y(T)  is  known  at  some  other  T.) 
Equations  (5),  (6a),  and  (8)  thus  also  serve  as  a  basis  for  the  discussion  of 
thermal  expansion  effects  in  electrical  transport.  Clearly,  if  K(T)  is  known 
one  can  use  Eqs.  (5)  and  (6)  without  further  approximation  to  describe  the 
thermal  expansion  part  of  p(T). 

Frequently,  the  pressure  coefficient  of  resistance  H  is  experimentally 

accessible  rather  than  3£np/3JlnV|  ,  One  then  employs 

T,V0 

II(T)  =  (3£nR/3P)  =  (  3 £nR/  3 &nV )  (3inV/3P) 

T,V0  T,V0  T,Vq 

=  -  (3JlnR/3£nV)  /BT(T)  (9a) 

TV 
1  }  v  o 

=  “  ((3£np/3JlnV)  -  l/3)/Bx(T)  (9b) 

T,Vn 
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or 


( 3£np/ 3£V)  =  -  BX(T)II(T)  +  1/3  (9c) 

TV  1 

L  > v  o 

where  R  is  the  resistance. 

We  will  refer  to  isochoric  and  thermal  expansion  parts  of  the  isobaric 
temperature  coefficients  of  resistivity  defined  as 

3£np  j 

r  = - 1  =  r  +  r  (ioa) 

3T  |t=8,P  vo  g 
where  the  isochoric  part  is  given  by 

T  =  (3R(T,V0)/3T)t=9  (10b) 

V0 

and  the  thermal  expansion  part  is  given  by 

T  =  (3Q(T)/3T)X=0  (10c) 

0 

=  0v(9)(3£np/32nV)  (lOd) 

Mo 

=  -  0V(  9)(bt(  e)rt(  0)  -  1/3)  (lOe) 

Note  that  the  description  of  thermal  expansion  effects  on  electrical 
transport  presented  here  is  model  independent.  That  is,  no  matter  which 
technique  is  used  to  compute  the  isochoric  relative  change  in  p  (i.e.,  R(T,V0) 
in  Eq .  (15)),  a  thermal  expansion  correction  Q(T)  is  required  to  obtain  the 
isobaric  p(T).  The  thermal  expansion  term  will  be  observable  in  p(T)  whenever 
II  (or  3£np/3£nV)  is  sufficiently  large  in  magnitude.  Equation  (6)  with  F(x) 
defined  in  Eq.  (4b)  obtains  if  V(T)  is  represented  in  terms  of  the  Gruneisen 
theory  of  thermal  expansion.  Generally,  K(T),  defined  in  Eq .  (6b),  is  a  weak 
function  of  T  in  comparison  to  F(T/0)  so  that  in  the  examples  given  here  we 
treat  K(T)  as  a  constant.  Clearly,  the  T-dependence  of  K(T)  or  a  more  general 
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theory  for  V(T)  and/or  ( 3£np/  3i.nV)  can  be  incorporated  into  the  formalism, 

T,V0 

in  a  straightforward  manner. 

RESULTS 

The  thermal  expansion  part  of  the  isobaric  p(T)  is  negligible  in 

comparison  to  the  isochoric  part  of  the  isobaric  p(T)  at  very  low  temperatures 

(T  <  9/15).  One  can  see  this  effect  in  Figures  2  through  5.  This  is  the  case 

because  F(T/6)  vanishes  as  (T/9)4  while  R(T,V0)  vanishes  approximately  as 

(T/0)2  at  very  low  temperatures.  The  thermal  expansion  term  is  approximately 

linear  in  T  above  the  Debye  temperature  and  yields  a  "correct ion"  to  the 

isochoric  temperature  coefficient  of  resistivity  T  whose  sign  is 

vo 

determined  by  the  sign  of  ( 3£np/ 3£nV) 

T,V0 

The  best  characterized  and  most  thoroughly  studied  low  resistivity 
(~  60  p^cra)  amorphous  alloy  is  a-MgyZng  (refs  4,7).  It  is  the  only  low 
resistivity  alloy  for  which  the  pressure  coefficient  of  resistance  H  has  been 
determined.  The  isobaric  p(T)  has  the  following  features:  (1)  a  small 
minimum  at  Tm  =»  7K;  (2)  a  maximum  at  Tj^  “  40K;  (3)  an  approximately  (T-T^)^/^ 
region  for  T  >  T^;  and  (4)  approximately  linear  T-dependence  near  300K  with  T 
“  -2 .0xl0_1+/K.  The  low  temperature  results  from  Reference  7  are  shown  in 
Figure  2. 

Since  a-MgyZng  is  the  only  low  resistivity  alloy  with  known  II  (ref  4)  ,  we 
will  provide  a  somewhat  extended  description  of  the  status  of  the  theoretical 
work  on  this  alloy  and  the  importance  of  thermal  expansion  in  determining 
p(T).  All  the  calculations  have  been  performed  in  the  context  of 
Bayra-Faber-Ziman  theory  (refs  8,9)  (the  diffraction  model). 
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Hafner  (ref  5)  computed  p(T)  in  a-MgyZn3  from  first  principles  in  a 
pseudopotential  based  procedure  with  no  saturation  corrections.  The  volume 
dependence  of  p,  H,  the  isobaric  T-dependence  of  p,  the  dynamical  structure 
factor,  and  the  pseudopotentials  were  all  computed  with  no  adjustable  or 
empirical  parameters.  The  agreement  with  the  data  is  remarkable.  Hafner 
concluded  that  thermal  expansion  effects  are  important  in  determining  p(T)  in 
(liquid  and)  amorphous  MgyZn3. 

Meisel  and  Cote  (ref  10)  described  isochoric  diffraction  model  studies  of 
p(T)  in  the  a-MgZn  alloys  (and  low  resistivity  alloys  in  general).  An 
effective  scattering  matrix  was  constructed  which  had  exclusively  s  and  p 
character,  satisfied  the  Friedel  sum  rule,  and  yielded  the  observed  magnitude 
of  p.  "Good"  quantitative  agreement  with  isobaric  data  was  obtained  when 
Pippard-Ziman  saturation  (refs  9,11-15)  was  included:  (1)  the  computed  T 
values  were  approximately  half  the  observed  values;  (2)  the  temperatures  T^ 
for  the  resistivity  maxima  were  overestimated  by  approximately  30  percent;  (3) 
the  sizes  of  the  maxima  ( p(T^j)/ p(0)-l)  were  overestimated  by  a  factor  of  two; 
(4)  the  temperatures  Tm  for  the  resistivity  minima  were  overestimated  by 
approximately  30  percent;  and  (5)  the  depths  of  the  minima  were  underestimated 
by  at  least  a  factor  of  two.  (These  results  are  closely  approximated  by  curve 
D  of  Figure  2.)  When  saturation  effects  were  not  included  (i.e.,  when 
standard  Baym-Faber-Ziman  theory  (ref  8)  was  applied)  qualitative  agreement 
with  the  data  was  still  obtained  (except  for  the  minimum  in  p) ,  although 
quantitative  agreement  with  the  data  was  severely  diminished.  (These  results 
are  closely  approximated  by  curve  A  of  Figure  2.) 
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A  later  discussion  of  isochoric  diffraction  model  transport  in  a-MgZn  by 
Cote  and  Meisel  (ref  14)  based  upon  phase  shifts  computed  for  Herman-Skillman 
(ref  15)  neutral  atom  wavef unctions ,  suggested  that  the  isobaric  p(T)  could  be 
fit  if  the  Debye  temperature  for  resistivity  was  taken  as  200K  instead  of  the 
thermal  value  of  295K.  (The  magnitude  of  p  computed  for  these  phase  shifts  is 
40  u^cm,  approximately  30  percent  less  than  the  measured  value.)  If  9  is 
taken  as  295K,  the  results  of  these  calculations,  shown  in  curve  D  of  Figure 
2,  are  not  significantly  different  from  those  of  the  "effective  potential" 
calculations  described  above. 

Figure  2  shows  results  of  diffraction  model  calculations  based  upon  the 
phase  shifts  employed  in  Reference  14.  Essentially  the  same  results  are 
obtained  in  constant  t-matrix  calculations  as  described  in  Reference  16  or 
effective  potential  calculations  as  described  in  Reference  10.  The  parameters 
used  In  the  calculations  are  listed  In  Table  I.  Isochoric  results  are  shown 
in  curves  A  and  D.  Isobaric  results  computed  using  the  pressure  coefficient 
of  resistivity  H  value  deduced  from  the  experimental  data  of  Fritsch  et  al 
(ref  4)  are  shown  in  curves  B  and  E  and  results  computed  using  the  H  value 
computed  by  Hafner  (ref  5)  are  shown  in  curves  C  and  F.  The  isochoric  part  of 
P  in  curves  A,  B,  and  C  is  the  standard  Baym-Faber-Ziman  (ref  8)  results, 
while  that  in  curves  D,  E,  and  F  include  Pippard-Ziman  saturation  (refs  9, 
11-13)  with  qj)A  =  17.1  where  is  the  Debye  wave  number  and  A  the  electron 
mean  free  path.  The  value  of  qpA  was  computed  in  the  free  electron 
approximation  (Eq.  (24)  of  Reference  16).  That  is, 

qDA  =  644(Z/2)1/3  (kpa^'^p/pSkn)'1  (11) 
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where  Z  is  the  electron  per  atom  ratio  and  ajj  the  Bohr  radius.  Although 
Figure  2  displays  results  for  0  <  T  <  105K,  our  discussion  pertains  to  results 
for  0  <  T  <  300K.  That  is,  when  we  refer  to  curve  F  we  mean  the  computed 
result  for  the  specified  parameter  values  for  0  <  T  <  300K. 

Curve  F  (incorporating  II  computed  by  Hafner  (ref  5)  and  qpA  =  17.1) 
agrees  with  the  data  quite  well:  (1)  Temperature  coefficient  of  resistivity  F 
is  given  within  10  percent;  (2)  the  position  of  the  maximum  p(T) ,  TM  is  given 
within  10  percent;  (3)  the  position  of  the  minimum  of  p(T),  Tm  is  given  within 
10  percent;  (4)  the  size  of  the  maximum  is  underestimated  by  about  30  percent; 
and  (5)  the  size  of  the  minimum  is  significantly  underestimated. 

Curve  E  (incorporating  H  deduced  from  the  data  of  Fritsch  et  al  (ref  4) 
and  q^A  =  17.1)  represents  almost  as  good  a  fit  to  the  data  as  curve  F:  (1)  F 
is  underestimated  by  about  25  percent;  (2)  T^  is  overestimated  by 
approximately  25  percent;  (3)  Tm  is  given  within  10  percent;  (4)  the  size  of 
the  maximum  is  given  within  10  percent;  and  (5)  the  size  of  the  minimum  is 
significantly  underestimated. 

The  curves  of  Figure  2  illustrate  the  action  of  thermal  expansion  effects 
for  positive  n.  They  also  illustrate  the  action  of  Pippard-Ziraan  saturation 
(refs  9,11-13)  for  a  low  resistivity  (qpA  ■=  17)  amorphous  metal.  Both  effects 
yield  more  negative  T,  shift  T^j  to  lower  values,  and  reduce  the  size  of  the 
maximum  of  p(T).  However,  at  very  low  T,  thermal  expansion  effects  become 
negligible  while  saturation  effects  are  still  important  (producing  minima  in 
p(T),  for  example).  (Note  that  Hafner' s  (ref  5)  transport  calculation,  which 
does  not  include  saturation,  produces  a  minimum  in  p(T) .  The  origin  of  the 
calculated  minimum  is  not  clear.  These  calculations  also  underestimate  the 
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size  of  the  minimum.)  It  is  evident  that  thermal  expansion  effects  are 
important  in  determining  p(T)  in  a-Mg7Zn3  and  are  of  the  same  order  as 
saturation  effects  for  this  alloy. 

When  thermal  expansion  effects  are  included  in  diffraction  model 
calculations  based  on  Heine-Aberenkov  pseudopotentials  for  Mg  and  Zn  tabulated 
in  Harrison  (ref  18),  qualitative  agreement  with  the  data  is  greatly  improved, 
but  quantitative  results  are  still  poor.  (However,  recall  that  Hafner's  (ref 
5)  pseudopotential  calculation  agrees  quite  well  with  the  data  for  a-MgyZng.) 

Figures  3  through  5  are  examples  of  the  effects  of  thermal  expansion  in 
determining  isobaric  p  vs.  T  curves  in  the  context  of  the  diffraction  model 
including  the  Pippard-Ziman  saturation.  The  parameters  and  variables  used  in 
the  figures  are  defined  as  follows: 

ci  =  12*(2kF2/MkB9  (12) 

where  M  is  the  averaged  ionic  mass  and  kF  is  the  Fermi  wavenumber.  The 
parameter  a  determines  the  scale  of  the  isochoric  variations  of  p  for  given 
values  of  2kp/kp,  qBA,  etc. 

For  practical  values  of  a,  R(T,V0)/a  is  approximately  independent  of  a. 
Thus,  we  define  the  normalized  relative  change  of  the  resistivity  as 


r  =  (p(T)/p0-l)/o  =  R(T,V0)/a  +  Q(T)/a  (13a) 

=  R(T,Vo)/a  +  w(T)F(T/9)  (13b) 

where,  in  the  approximation  scheme  which  treats  K(T)  as  constant, 

w(T)  «  w  =  w(  9)  =  -1.05  9Sv(9)n/a  (13c) 

The  thermal  expansion  part  of  T  can  then  be  expressed  as 

T  =  0.95  a  w(9)/9  (13d) 
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We  also  define  a  normalized  temperature, 

t  =  T/0  ( 13e) 

Thus,  Figures  3  through  5  are  plots  of  normalized  relative  change  of  the 
resistivity  r  versus  normalized  temperature  t  for  specific  values  of  2kp/kp 
and  q^A  and  a  range  of  values  for  the  parameter  w.  The  w  values  (viz.,  0, 

±  0.2,  and  ±  0.5)  were  chosen  to  be  representative  of  the  range  to  be  found  in 
amorphous  metals.  The  curves  can  be  linearly  interpolated  or  extrapolated  to 
obtain  results  for  other  w  values  or  to  allow  for  a  T-dependent  w(T).  For 
example,  in  a-MgyZng,  w  a  0.16  for  the  experimental^  II  and  0.33  for  the 
theoretical  (ref  5)  H. 

Results  based  upon  potassium  pseudoatom  phase  shifts  computed  by  Young, 
Meyer,  and  Kilby  (ref  19)  for  2kp/kp  =1.1  and  qgA  =  12  are  shown  in  Figure  3. 
The  pseudoatom  phase  shift  absolute  squared  t-matrices  are  very  similar  to 
absolute  squared  pseudopotential  matrix  elements  except  that  they  do  not 
exhibit  perfect  cancellation  near  2kp.  Thus,  Figure  3  is  very  similar  to  a 
pseudopotential  result.  Isochoric  p(T,V0)  has  a  maximum  near  9/2  and  r  is 

Vq 

slightly  negative.  Thermal  expansion  effects  of  the  order  of  magnitude  of 
that  measured  by  Fritsch  et  al  (ref  4)  yield  isobaric  p(T)  which  are  in 
reasonable  agreement  with  the  forms  observed  in  simple  metal  amorphous  alloys. 
Note,  in  particular,  the  dramatic  reductions  in  the  size  of  the  maximum  in 
p(T)  and  its  movement  to  lower  T^  as  w  varies  from  0  to  -0.5. 

Figure  4  shows  results  computed  in  the  constant  t-matrix  approximation  as 
described  in  Reference  16,  at  2kp  =  1.15  kp  and  qpjA  =  12.  The  isochoric  and 
negative  w  curves  in  this  figure  are  not  significantly  different  from  those 
computed  for  the  Mg  and  Zn  phase  shifts  with  2kp  =  1.11  kp,  and  qpA  =  17.1 
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shown  in  Figure  2.  This  figure  illustrates  the  reduction  in  magnitude  and 
motion  of  the  maximum  in  p(T)  as  w  varies  from  +0.2  to  -0.5.  For  w  values 
slightly  more  negative  than  -0.5,  p(T)  is  monotonic  decreasing  with  increasing 
T  and  for  w  =  +0.5,  the  maximum  in  p(T)  is  gone  and  F  is  positive. 

Figure  5  shows  results  in  the  constant  t-matrix  approximation  (as  in 
Figure  4)  with  2kp  =  1.3  kp  and  qpjA  =  6.  These  curves  are  thus  appropriate  to 
larger  electron  per  atom  ratio  z  and  resistivity  than  those  shown  in  Figures  2 
through  4.  We  would  not  consider  a  calculation  for  q^A  =  6  as  representative 
of  low  p  alloys.  The  isochoric  p(T)  has  a  small  negative  T  and  a  maximum  near 
0/4,  the  positive  w  curves  have  positive  T  and  no  maxima,  the  curve  for  w  = 
-0.2  has  a  small  maximum  near  0/6  and  negative  T,  and  the ' curve  for  w  =  -0.5 
is  monotonic  decreasing. 

The  curves  shown  in  Figures  3  through  5  are  not  intended  to  provide  an 
exhaustive  survey  of  the  possibilities  for  the  inclusion  of  thermal  expansion 
effects  in  p(T)  by  any  means.  For  example,  each  of  the  sixty  isochoric 
results  presented  in  Reference  16  could  be  the  basis  for  a  set  of  isobaric 
results  for  different  w(T).  Thus,  the  present  selection  is  necessarily 
somewhat  arbitrary. 

An  interesting  case,  which  we  have  not  shown  in  the  figures  is  p  versus  T 
for  extreme  saturation  (qpA  >  1) .  It  is  possible  to  select  positive  w  values 
in  this  limit  which  yield  p  versus  T  curves  which  resemble  those  seen  in  very 
high  resistivity  (p  >  200  pWcm)  metals  (ref  20).  However,  the  choice  of 
parameters  (a  and  w)  to  fit  observed  p(T)  seems  to  be  too  limited  to  explain 
an  apparently  universal  form.  Perhaps  realistic  T-dependent  w(T)  could 
explain  the  extreme  saturation  p(T)  forms. 
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Another  interesting  case  not  illustrated  is  the  moderate  or  low  p  alloy, 
whose  isochoric  p(T)  has  T  >  0.  For  such  alloys  (with  positive  T)  inclusion 
of  thermal  expansion  can  yield  the  typically  observed  negative  T  with 
appropriate  maximum,  etc. 

Let  us  now  address  the  general  question  of  thermal  expansion  effects  in 
actual  amorphous  alloys.  Let  TM  stand  for  a  transition  metal  constituent,  M 
stand  for  a  non-TM  metallic  constituent,  and  NM  a  non-metallic  constituent. 
Fritsch  et  al  (ref  4)  made  the  following  observations  concerning  the  pressure 
coefficient  of  resistivity  T  in  amorphous  alloys:  (1)  M-TM  alloys  exhibit 
negative  H;  (2)  TM-TM  alloys  exhibit  negative  H;  and  (3)  NM-TM  systems  exhibit 
negative  or  small  H.  Essentially  the  same  conclusions  were  drawn  by  Cote  and 
Meisel  (ref  2)  (based  on  considerably  less  data)  in  their  earlier  study  of 
pressure  effects  in  electrical  transport.  However,  in  contrast  to  Cote  and 
Meisel  (ref  2),  Fritsch,  et  al  (ref  4)  also  concluded  that  there  is  no 
correlation  between  the  H  and  the  magnitude  of  p.  It  appears  to  us  that  the 
data  of  Fritsch,  et  al  (ref  4)  are,  in  fact,  consistent  with  our  conclusion 
that  the  envelope  of  H  versus  p  values  becomes  narrower  at  higher  p  values, 
although  the  narrowing  is  not  as  rapid  as  estimated  in  Reference  2.  Also,  one 
should  note  that  the  higher  p  amorphous  alloys  tend  to  be  less  compressible 
than  the  low  p  alloys  so  that  the  narrowing  of  the  envelope  of  the  II  versus  p 
data  may  be  equivalent  to  a  p  independent  envelope  for  ( 3 in p/ 3 &nV ) 

T  V* 

L  > v  o 

versus  p. 

No  general  conclusions  concerning  II  in  low  p  or  M-M  amorphous  alloys  can 
be  drawn  from  extant  pressure  data.  However,  there  are  reasons  to  believe 
that  H  may  be  positive  or  small  for  M-M  amorphous  alloys:  (1)  the  only  system 
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studied  (a-Mg7Zn3)  exhibits  a  relatively  large  positive  II;  (2)  values  of  qoA 
smaller  than  those  indicated  by  Eq .  (11)  and/or  significantly  smaller  9  than 
the  thermal  values  are  required  in  the  context  of  the  diffraction  model  to 
yield  isochoric  p(T,VQ)  which  are  in  agreement  with  the  extensive  body  of 
isobaric  p(T)  data  in  M-M  amorphous  alloys  determined  by  Mizutani  and 
coworkers  (refs  7,21).  The  discrepancies  between  the  observed  isobaric  p(T) 
and  the  computed  isochoric  p(T,V0),  without  adjusted  qpA  and  6  are  consistent 
with  a  thermal  expansion  contribution  appropriate  for  positive  H.  (This 
picture  appears  to  hold  together  in  detail  for  a-MgyZng.) 

Thus,  the  curves  in  Figures  3  through  5  for  w  >  0  are  representative  of 
high  p  (p  >  100  pQcm)  alloys,  which  are  generally  TM-TM,  M-TM,  or  NM-TM 
systems;  and  the  curves  in  Figures  2  through  5  for  w  <  0  are  (tentatively) 
representative  of  low  p  (p  <  100  yOcm)  alloys,  which  are  generally  M-M 
systems.  We  keep  in  mind  that  in  general  w(T)  will  be  T-dependent  and  allow 
for  the  possibility  of  T-dependent  linear  interpolation  or  extrapolation  of 
these  curves.  The  form  of  F(T/9)  shown  in  Figure  1  can,  of  course,  be 
combined  with  arbitrary  isochoric  p(T,VQ)  (or  R(T,V0)/a)  and  appropriate  K(T) 
(or  w(T) )  to  obtain  isobaric  p(T)  (or  r(T))  curves. 

CONCLUSIONS 

1.  Thermal  expansion  effects  are  important  in  determining  isobaric  p(T) 
in  amorphous  metals.  Such  effects  are  independent  of  the  particular  treatment 
appropriate  to  the  isochoric  problem.  For  example,  thermal  expansion 
contributes  at  least  25  percent  of  the  observed  isobaric  temperature 
coefficient  of  resistivity  T  in  a-MgyZng. 
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2.  All  TM-TM,  M-TM,  and  NM-T,  amorphous  alloys  studied  (which  generally 
have  p  >  100  pftcm)  have  small  or  negative  pressure  coefficient  of  resistivity 
n  and  hence,  from  Eq.  (13),  have  small  or  positive  w(T). 

Values  of  II  are  not  available  for  M-M  and/or  low  resistivity  amorphous 
metals,  with  the  notable  exception  of  a-MgyZng.  However,  the  preponderance  of 
relatively  large  negative  isobaric  T  values  in  the  alloys  studied  (refs  7,21) 
and  the  observation  in  a-MgyZng  suggests  that  positive  or  small  H,  and  hence 
negative  or  small  w(T) ,  may  be  appropriate  for  these  low  resistivity  alloys. 

3.  Alloys  having  negative  T  exhibit  either  monotonic  decreasing  p(T) 
or  a  maximum  in  p(T).  The  p(T)  exhibiting  a  maximum  often  exhibit  a  minimum 
in  p(T)  at  a  lower  T  as  well.  In  some  of  the  cases  where  a  minimum  is  not 
observed,  the  possibility  of  the  observation  of  a  minimum  is  precluded  by  lack 
of  experimental  precision  or  a  transition  to  the  superconducting  state. 

Monotonic  decreasing  p(T)  is  usually  observed  in  high  p  cases  (refs  9,20) 
and  the  maximum  (with  or  without  minimum)  in  p(T)  is  usually  observed  in  low  p 
or  moderate  p  cases.  (We  suggest  that  these  forms  for  p(T)  result  from 
Pippard-Ziman  saturation  (refs  9,11-13)  with  appropriate  q^A.) 

The  effect  of  thermal  expansion  on  P(T)  having  a  maximum  for  positive 
(negative)  w(T)  is  to  increase  (decrease)  the  size  of  the  maximum  and  to  move 
the  maximum  to  higher  (lower)  temperature  than  given  by  isochoric  theory.  It 
is  also  possible  for  thermal  expansion  with  positive  w(T)  to  produce  isobaric 
p(T)  with  minimum  and  maximum  for  monotonic  decreasing  isochoric  p(T,VQ)  and 
vice  versa.  (Note  Chat  thermal  expansion  described  by  Gruneisen  theory,  i.e., 
Eqs.  (3)  and  (4) ,  has  negligible  effect  at  the  temperatures  where  minima  have 
been  observed.) 
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Thermal  expansion  effects  with  positive  w(T)  can  produce  upward  curvature 
for  T  >  9/5  in  monotonic  decreasing  p(T)  curves.  Such  curvature  is  observed 
in  some  high  p  alloys. 

4.  Positive  T  alloys  are  predominantly  of  low  resistivity  (p  <  100  uffcm) 

and  one  can  expect  to  find  upward  (downward)  curvature  in  p(T)  for  T  >  0/5  for 

w(T)  positive  (negative).  We  have  also  shown  that  thermal  expansion  effects 

for  positive  w(T)  can  yield  isobaric  p(T)  with  positive  T  for  isochoric 

P(T,V0)  having  maximum,  minimum,  and  negative  F  and  vice  versa. 

vo 

5.  Thermal  expansion  effects  are  important  in  determining  p(T)  in  the 

specific  case  of  a-Mg7Zn3-  Good  agreement  with  the  data  is  obtained  with 

Gruneisen  thermal  expansion  theory,  Eqs.  (3)  and  (4),  with  w(T)  assumed 
constant  equal  to  w(9)  empirical  or  thermal  values  of  parameters  and 
diffraction  model  transport  incorporating  Pippard-Ziman  saturation  with  q^A  = 
17.1  and  scattering  phase  shifts  based  on  Herman- Ski liman  neutral  atom 
wavef unctions  (ref  15).  Hafner  (ref  5)  also  obtained  good  agreement  with  the 
a-MgyZn3  resistivity  data  in  a  first  principles  pseudopotential  based  standard 
(i.e.,  no  saturation)  diffraction  model  calculation  which  incorporated  thermal 
expansion  effects. 

Besides  serving  to  demonstrate  the  importance  of  thermal  expansion 
effects  in  electrical  transport  in  amorphous  metals,  the  a-Mg7Zn3  results 
raise  some  disturbing  questions  concerning  diffraction  model  procedures.  For 
example,  the  low  temperature  minimum  in  p(T) ,  which  in  our  calculations  is  a 
saturation  effect,  is  a  consequence  of  some  feature  of  the  standard 
diffraction  model  in  Hafner' s  work.  We  do  not  know  which  features  of  Hafner' s 
work  are  significantly  different  from  ours.  Is  Hafner' s  self-consistent 
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phonon  spectrum  significantly  different  from  the  Debye  form?  Is  his  treatment 
of  thermal  expansion  significantly  different  from  the  simple  Gruneisen 
approach  adopted  here?  These  (and  other)  questions  are  presently  unresolved. 
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TABLE  I.  PARAMETERS  EMPLOYED  IN  THE  a-MgyZn3  CALCULATIONS 


0  =  295K  is  from  Mizutani  and  Mizoguchi  in  Reference  7. 

2kp/kp  =  I. II  where  kp  is  from  Matsuda  and  Mizutani  in  Reference  7. 

and  kp  is  from  Mizoguchi  et  al  in  Reference  7. 

9D  =  anc*  single  branch  Debye  spectrum  are  assumed. 
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|  89.1 
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103  II  (kbar-1) 

1 

|  14.00 
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i 

i 

!  9.17 

12.6 

1 

1 

!  2.39 

1 

4.8 

*These  values,  which  are  computed  by  the  "law  of  mixtures",  and  the  pressure 
coefficient  of  resistivity  II  values  shown  are  used  in  the  transport 
calculations . 
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Figure  2.  The  relative  change  in  resistivity  versus  temperature  for  a-MgyZn^  based  on  the 

parameters,  in  Table  I.  Curves  A,  B,  and  C  are  computed  for  q^A  =  <»  (the  standard 
diffraction  model);  curves  D,  E,  F  for  q^A  =  17.1.  Curves  A  and  D  are  isochoric 
results.  Curves  B  and  E  include  thermal  expansion  effects  with  the  measured  (ref  4) 
CpCR-  Curves  C  and  F  include  thermal  expansion  effects  with  Ilafner's  (ref  5) 
theoretical  Cp£R. 


0.04- 


Figure  3.  Normalized  resistivity  difference  versus  normalized  temperature  for  potassium 
pseudoatom  phase  shifts  (ref  19)  with  ?.kF/kp  =  1.1  and  qDA  =  12.  The  p 
values  are  0.5,  0.2,  0,  -0.2,  and  -0.5  for  curves  A,  B,  C,  D,  and  E,  respectively. 
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